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Abstract 

We study the production of prompt diphotons in the central region of rapidity within the frame- 
work of the quasi-multi-Regge-kinematics approach applying the hypothesis of quark and gluon 
Reggeization. We describe accurately and without free parameters the experimental data which 
were obtained by the CDF Collaboration at the Tevatron Collider. It is shown that the main con- 
tribution to studied process is given by the direct fusion of two Reggeized gluons into a photon pair, 
which is described by the effective Reggeon-Reggeon to particle-particle vertex. The contribution 
from the annihilation of Reggeized quark-antiquark pair into a diphoton is also considered. At 
the stage of numerical calculations we use the Kimber-Martin-Ryskin prescription for unintegrated 
quark and gluon distribution functions, with the Martin-Roberts-Stirling-Thorne collinear parton 
densities for a proton as input. 
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I. INTRODUCTION 



The study of production of photons with large transverse momenta in the hard inter- 
action between two partons in the high-energy hadron collisions, so-called prompt photon 
production, provides precision tests of perturbative quantum chromodynamics (QCD) as 
well as information on the parton densities within a hadron. Also, these studies are our 
potential for the observation of a new dynamical regime, namely the high-energy Regge 
limit, which is characterized by the following condition v^S* >> /x >> Aqcd, where ^/S is 
the total collision energy in the center of mass reference frame, Aqcd is the asymptotic scale 
parameter of QCD, /x is the typical energy scale of the hard interaction. At this high-energy 
limit, the contribution from the partonic subprocesses involving t— channel parton (quark or 
gluon) exchanges to the production cross section can become dominant. In the region under 
consideration, the transverse momenta of the incoming partons and their off-shell properties 
can no longer be neglected, and we deal with "Reggeized" t— channel partons. In this kine- 
matics, so-called quasi-multi- Regge kinematics, the particles (multi-Regge) or the groups of 
particles (quasi-multi- Regge) produced in the collision are strongly separated in rapidity. In 
our discussion, the both photons from the photon pair are produced in the central region of 
rapidity. 

n 

The quasi-multi-Regge-kinematics (QMRK) approach [1[ is particularly appropriate for 
this kind of high-energy phenomenology. It is based on an effective quantum field theory im- 
plemented with the non-Abelian gauge-invariant action including fields of Reggeized gluons 
2] and Reggeized quarks In the QMRK approach we can do calculations with Reggeized 
quarks l^Q] that presents an open question in the other uncollinear factorization scheme, 
namely kx— factorization approach 7|, following which we can operate correctly with off- 
shell gluons only. Recently it was shown also, that the calculation in the next to leading 
order (NLO) in strong coupling constant within the framework of QMRK approach can be 
done [8]. 

Our previous studies of charmonium and bottomonium production Q], open charm pro- 
duction [5|and inclusive prompt photon production at the DESY HERA and at the Fermilab 



Tevatron [6| demonstrated the advantages of the high-energy factorization scheme based on 
the QMRK approach over the collinear parton model as far as the description of experimental 
data is concerned. 
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This paper continues our study of the inclusive prompt photon production at Tevatron 



which was performed recently 



6| in the framework of the the QMRK approach. We consider 



here the diphoton production at the Fermilab Tevatron in the central region of rapidity 



applying the hypothesis of quark and gluon Reggeization 

It is shown that the main contribution to studied process is given by the direct fusion 
of two Reggeized gluons into a photon pair, which is described by the effective Reggeon- 
Reggeon to particle-particle vertex. The contribution from the annihilation of Reggeized 
quark - antiquark pair into a diphoton is not small and it is taken into account additionaly. 
We do not take into consideration the contribution of the fragmentation mechanism of the 
prompt photon production which is strongly suppressed by the isolation cone condition 

At the stage of numerical calculations we use the Kimber-Martin-Ryskin prescription 
for unintegrated quark and gluon distribution functions, with the Martin-Roberts-Stirling- 
Thorne (MRST) collinear parton densities for a proton as input . 

This paper is organized as follows. In Sec. Dll the relevant Reggeon-Reggeon to particle- 
particle effective vertices are presented and discussed. In Sec. IHIl we describe diphoton 
production spectra at the Tevatron Collider. In Sec. HV] we summarize our conclusions. 

II. BASIC FORMALISM 

In the phenomenology of strong interactions at high energies, it is necessary to de- 
scribe the QCD evolution of the parton distribution functions of colliding particles start- 
ing with some scale fiQ, which controls a non-perturbative regime, to the typical scale /i 
of the hard-scattering processes, which is typically of the order of the transverse mass 



Mt = y M 2 + |/ct| 2 of the produced particle with mass M and transverse momentum kr- In 
the region of very high energies, in so-called Regge limit, the typical ratio x = /i/ becomes 
very small, s < 1. That leads to large logarithmic contributions of the type [a s ln(l/x)] n 
in the resummation procedure, which is described by the Balitsky-Fadin-Kuraev-Lipatov 



(BFKL) evolution equation 12] or other BFKL-like ones for unintegrated gluon (quark) 
distribution functions |q T | 2 , /i 2 ). Correspondingly, in the QMRK approach [1], the 

initial-state t-channel gluons and quarks are considered as Reggeons, or Reggeized gluons (R) 
and Reggeized quarks (Q) . They are off-mass shell and carry finite transverse two- momenta 
qr with respect to the hadron beam from which they stem. 
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The advantages of the QMRK approach include: first, it uses gauge-invariant amplitudes 
and is based on a factorization hypothesis that is proven in the leading logarithmic approx- 
imation; second, it carries over to non-leading orders in the strong-coupling constant, as 
recently proven third, it works both with Reggeized gluons and with Reggeized quarks. 
The Reggeization of amplitudes provides the opportunity to efficiently take into account 
large radiative corrections to processes in the Regge limit beyond what is included in the 
collinear approximation, which is of great practical importance. The particle Reg geiz ation 
is known effect in the high-energy quantum electrodynamics (QED) for electrons [13|, and 
for gluons [12J and quarks 3|, |9[] in QCD. 

Recently, in Refs. [3^, Il4 1 , the Feynman rules for the effective theory based on the non- 
Abelian gauge— action includi „ g ne.ds of Re gg e ized glu0nS Q a nd Re ggelZ ed qUarks 
[9I were derived for the induced and some important effective vertices. The effective vertices 
for the 2-^2 processes with Reggeized gluons in the initial state only, C ~* 99 and C "* qq , 



were obtained in Refs. jl2|, |l5] , the effective vertices C QQ ^ 99 and C R ^ 9q with Reggeized 
quark and Reggeized gluon in the initial state were obtained in Ref. [sj]. For our purposes we 
need to construct the effective vertices with two photons in the final state, i.e. C^^ 77 and 
(jRR^-yy _ rp^g g rs ^ one can De eaS y obtained from the above mentioned vertices after cutting 
the contribution from the three-gluon vertex and replacement the quark-gluon vertex by the 
same quark-photon vertex. The effective vertex C^ - * 77 can be obtained using the Feynman 



rules 14j for the Reggeized gluons and the well-known fourth-rank vacuum polarization 
tensor, which corresponds to the set of quark-box diagrams [l6]. 

We perform calculations in the laboratory frame of the Tevatron Collider, where proton 
and antiproton beams have equal energies, E p = E p = 900(980) GeV. We introduce two light- 
cone vectors corresponding to the massless particle momenta as follows: Pi = Pi (1,0,0, 1) 
and P 2 = P 2 ( 1,0,0 — 1), where Pi = E p and P 2 = Ep. Also we define two additional 
four- vectors (n + ) At = P^ /P 2 and (n - ) M = P± /E\. For any arbitrary four-momentum k^, we 



define k^ 1 = k ■ n ± = k^n^. We use the vertex functions defined in Ref. [3|, which describe 
transition of the Reggeized quark with the four-momentum q and gluon or photon with the 
four-momentum k into the on-shell massless quark with four-momentum k + q: 



7i +) (A:,g)=7, + g^ = 7, + g^, (i) 

n7, A P 



V(k,q)= 1 , + q^= 1 , + q 1 ^. (2) 



The effective Reggeon-Reggeon to particle-particle vertex C®Q~* ry (qi, q 2 , ki, k 2 ), which 
describes the Reggeized quark - Reggeized antiquark annihilation into a photon pair in the 
process 

Q(qi)+Q(q2)^ 7(fci) +J{h), (3) 

can be presented as follows: 



(?2 + &2, 



+ 7< +) (*i, g 2 ) ( ^ + gi) + A^fa, -ft)] , (4) 

where momenta of the initial Reggeized parton are qn )2 ) = £(i,2)-P(i,2) + <?(i,2)t, an d the 
induced term A^(gi, g 2 ) has the form 

A^fai, gjj) = gi-pr- + 92ttt+ • (5) 
The vertex (jlj) satisfies the gauge-invariant condition, which reads as follows 

C% Q ^(qi, ft, k h k 2 )k% = C%^(q u q 2 , k h k 2 )k\ = 0. (6) 

The effective gauge-invariant vertex C RR ~*' y ' y (qi,q 2 ,k 1 ,k 2 ), which describes the direct 
Reggeized gluon fusion into a diphoton in the process 

R(qi) + R(q 2 )^l(h) +1 (k 2 ), (7) 

can be presented as follows 

C R ^{ qi) q 2 MM) = U { - )a (q 1 )U { + )l3 (q 2 )G a ^ I/ (q 1 ,q 2 ,k 1 ,k 2 ), (8) 

Gap^qi, q 2 ,ki,k 2 ) is the vacuum polarization tensor of forth rank corresponding the set of 
quark-box diagrams, which are taken with the relevant color factor. 

It is suitable to define the new csalar vertex C 77 (q±, q 2 ,ki, k 2 ), in which the summation 
over final photon polarizations has been performed: 

C RR ^(q 1 , q 2 , h, k 2 ) = £ e^h, M)e v {k 2 , \ 2 )C™-*n (q u q 2 , h, k 2 ), (9) 

Ai,A 2 

where Ai) and e u (k 2 , X 2 ) are the polarization four-vectors of final photons. The vertex 

© is presented as follows 

C RR ^{q x , q 2 , h, k 2 ) = n ( T )a (qi)n^(q 2 )G a0 ( qi , q 2 , h, k 2 ), (10) 



where 



G a/3 (q x ,q 2 ,k x ,k 2 ) = £ft (h, X x )e u (k 2 , X 2 )G a p flv (qi, q 2 , h, k 2 ) 



The exact expression for the tensor G a p{q X: q 2 , k x , k 2 ) have been obtained in the Ref. [16 ]. 
We use massless four-quark scheme to calculate this tensor. 

To obtain the vertex (TlTJT) we operate with projectors on the Reggeized gluon states which 
can be presented in the two equivalent forms U^ u (q 2 ) = ^ 2T or U^ u (q 2 ) — ' — 



\Q2T\ \Q2T\ 

and H^ u (qi) = 1 1T or H^ u (q x ) = 1 ^ j— . On the contrary to the definition used 

|9it| ]<?it| 



in Ref. [3j, we do not include on-shell quark spinors in our equations for the effective vertex 
OH). We also used different normalization for the projector IT^ , which reads accordingly 
Ref. as = (n^^. Our definition implies that the squared Reggeized amplitude 

in the QMRK approach is normalized to the squared parton amplitude for on-shell quarks 
and gluons when q XT = q 2T = 0. 

At the next step we write the squared matrix elements of the above mentioned Reggeized 
parton processes taking into account kinematical conditions of the Tevaron Collider. 

The squared matrix element for the direct diphoton production via annihilation ([3]) of 
Reggeized quark from a proton (gi = x x P x + q X T) and Reggeized antiquark from a antiproton 
(q 2 = x 2 P 2 + q 2 x) is obtained from the effective vertex (jlj) and it is presented as follows 

|M(Q P Q P -^77)| 2 = f^a 2 + w x S + w 2 S 2 + w 3 S 3 ), (11) 

o a x a 2 b x b 2 btu v ' 

where a x = 2k x ■ P 2 /S, a 2 = 2k 2 ■ P 2 /S, b x = 2k x ■ PJS, b 2 = 2k 2 ■ P x /S, S = 2P X ■ P 2 , 
t\ = -qf T , t 2 = —q 2T , s = (q x + q 2 ) 2 , t = (q x — k x ) 2 , u = (q x — k 2 ) 2 , and we apply that 
s + i + u = —t x - t 2 , x x — a x + a 2 , x 2 = b x + b 2 , 

w = t 1 t 2 (t x +t 2 ) ~iu(i + u), (12) 

-w x = t x t 2 (a x - a 2 )(b x -b 2 ) + t 2 x x (b 2 i + b x u) + 

+ t x x 2 (a x i + a 2 u) + iu(a x b x + 2a 2 b x + 2a x b 2 + a 2 b 2 ), (13) 

— w 2 = b x b 2 x\t 2 + a x a 2 x\t x + a x b 2 t(x x b x + a 2 b 2 ) + a 2 b x u(a x b x + a 2 x 2 ), (14) 
-w 3 = a x a 2 b x b 2 {a x b 2 [^j + a 2 &i(j)). (15) 
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When we consider the collinear approximation, in which one has ai 



u 



x 2 S 



a 2 



x 2 S 



t 



u 



h X 1 X 2 



X\b Xib o 

on-shell parton amplitude is obtained 



, and ti = t 2 = 0, the well-known answer for the squared 



|M(gg->7 7 )| J 



32 2 4 2 ft U\ 
3 H x u t' 



(16) 



The squared matrix element for the symmetric subprocess Q p Qp — > 77 is treated very 
similarly, one can be obtained from f[TT|) by the replacements: ti <->• t 2 , t «-> u, ai <-» 61, 
a 2 62, arid 07 <-» It is easy to see that these replacements do not change formula (111ft . 



The squared matrix element \M(RR — > 77) | 2 for the direct diphoton production via the 
fusion of Reggeized gluons ([7j) is obtained from the effective vertex (1101) . The analytical 
answer for the squared matrix element is known, see Ref. l?j], in the following approximation: 
m q = and ti >2 <C s, t, u, where m q is the mass of an quark in the relevant amplitude. It 



can be presented in the following form: 



\M{RR^^)\ 2 = 2a 2 a 2 (^e 2 q ) (h + f 2 + h) 



l f x 2 
l 1 



i=l 



(17) 



where n/ = 4 is the number of active quark flavors, 



fi 



's 2 + t 2 



+ 



+ 



's 2 + u 2 



log^ 



— s 

t 



2—\og( 

u \ t 



+ 



t 2 



,, -s\ s — u, f—s 
log" — +2^- log — 
u J t \ u 



+ 



u 



h 



+ 



s 2 + i 2 



7T" 



2— log(-) 



s 2 + u 2 



>og'[-) + 



log 



i 2 + u 2 



log 



u 



+ 



log 2 



7T 



-T^logf^ 



f - n 

h - y 



s 2 + g 

u 2 



+ 



s 2 + u 2 



1os It) + £ t 



+ 4. 



The formula (JT7J) formally coincides with the result, which is obtained in the collinear par- 
ton model. Any way, we incorporate the (fTTI) with the off-shell kinematics for the initial 
Reggeized-gluons. Of course, the exact squared matrix elements with the off-shell initial 
particles is needed. Our relevant calculations for arbitrary values of t\ and t 2 are in the 
progress. 



III. PROMPT DIPHOTON PRODUCTION AT TEVATRON 



During last decade the CDF and DO Collaborations at Tevatron obtained experimental 
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19 1, in association with a jet 



221 ] . The inclusive prompt 



data for the production of prompt photons: inclusively 
[20I . in association with a heavy quark (c, 6) {2]]], and in pair 
photon production was studied in the QMRK approach in the Ref. It was shown that 
the main mechanism of the inclusive prompt photon production in pp collisions is the fusion 
of Reggeized quark and Reggeized antiquark into a photon, via the effective Regeon-Reggeon 
to particle (to photon) vertex C®®^ 1 '. We have described well the inclusive photon transverse 
momentum spectra measured by the CDF and DO Collaborations [181 . 1191 ] in the wide region 
of the photon transverse momentum and pseudorapidity. 

Here we consider the production of prompt photon pair or diphoton in the framework 
of the QMRK approach and compare our predictions with the CDF data [22J. The cross 
sections are measured differentially as a function of the diphoton transverse momentum 
(Pt), the diphoton invariant mass (M), and the azimuthal angle between the two photons 
(A<£>). In this experimental analysis, the isolation condition required that the transverse 
energy sum in a cone of radius R = 0.4 (in ip — r\ space) about the photon direction, minus 
the photon energy, be less than 1 GeV. It is important to note, that the both photons in 



the CDF experiment 



22] are produced in the central region of rapidity (pseudorapidity) 



1 77^2 1 < 0.9. We can consider such photons as a quasi-multi-Regge group of particles. 

There are two mechanisms of prompt photon production: the production of direct pho- 
tons and the production of photons in fragmentation of produced quarks and gluons into 
photons. We have obtained [6y that the contribution of fragmentation mechanism is strongly 
suppressed by the isolation cone condition, which is applied to the experimental data for 
the inclusive photon production Q]- We estimate one as small, about 5 %, and do not 
take into account this contribution in the presented analysis too. 

The leading contributions in diphoton production are originating from the 2^2 pro- 
cesses ((3]) and ((71). The first one is order of a 2 , the second one is order of a 2 a 2 s . 

The factorization formulas for the hadron cross sections in the QMRK approach are the 
following: 

MPP - T,X) = ^L MQ Q 77 ) x 

J X\ J 7T J X 2 J 7T 



s 



X 



®q/p(xi, h, /i 2 )$Q/p(x 2 , t 2 , fi 2 ) + $Q/ p (Xi, h, fi 2 )$ Q/p (x 2} t 2 , /x 2 ))l8) 
da{vv - llX) = [*Ei[*M[**L[ d ^L d a(RR - 77 ) x 

J X\ J TT J X 2 J 7T 

x$ fl /p(xi, ti, fi 2 )^ R/p (x 2 , t 2 , ,u 2 ) (19) 

where $q jj/ p (x, t, /i 2 ) are the unintegrated parton distribution functions, da(RR, QQ — > 77) 
are the Reggeized parton cross sections. The unintegrated distribution functions and the cor- 
responding collinear F gtg / p (x, fi 2 ) distribution functions are connected by the normalization 
condition 

xF g , g/p (x, n 2 ) = / ®Q !R / P (x,t,fi 2 )dt, (20) 

•J 

that ensures the correct transition to the collinear-parton limit of equations flTSJ) and ffT9l . 
The Reggeized parton cross sections are defined in the following manner 



da(RR, QQ -> 77) = (2vr) j (g x + q 2 - fa - k 2 ) [[ . , (21) 

I RR,QQ i=\ { A7r > LK i 

where Irr^qq = 2x\X 2 S is the flux factor of the incoming particles, kl = (k®, fax, kf) and 
^2 = {k 2 , k-xr-, k 2 ) are the four-momenta of produced photons. 

Using the formulas ffTS]) . (fT9~l) and f[2"Tj) . we obtain the differential cross sections for the 
diphoton production, which can be written in case of the process (J7J) as follows: 



daipp^^X) 1 /• /■ /•„ /• /■ /• fc 1T fc 2r |M(i2i2-f 7 7)| 2 



<iti y d<^i y d^ir y y d&2T y ^2 



d(Ayj) 16vr 3 

x®r/p(xi, h, /i 2 )$R/ P -(x 2 , t 2 , /i 2 ) (22) 
dcr(pp -> 77X) p T 



dpT 167T 3 



<iti y d(px J dfax J dr/i J dk 2 r J dr] 2 



x^ R/p (x u t u /i 2 )$ i?/p -(x 2 , t 2 , /j 2 ) (23) 



'r/pY-zi, <>i,V> )VR/p 
da{pp — > 77X) M 



c?ti y dyji y d/ciy y <ir/i y dk 2 T J drj 2 



(x 


ix 2 S) 2 


\M(RR 


_> 77 )|2 


(xix 2 S) 2 \ 


sin(A(^) 


\M(RR 


_> 77 )|2 


(x l x 2 S) 2 \ 


sin(A<^) 



dM 16tt 3 

X$fl/ p (x ls ti, /i 2 )$R/ P -(x 2 , £ 2 , yu 2 ) (24) 

where fci j2 T = |^i,2t|, 771,2 are the photon pseudorapidities, <p\ is the azimuthal angle between 
£ 1T and<f ir , 2* = (k°+k°+kf+kl)/VS, x 2 = (fc?+fcjj-A£-A2)/VS, k% = ^(e^+e"^ 2 ), 
k\ 2 = -^p-^ 1 ' 2 — e _m ' 2 ). The differential cross sections for the diphoton production in case 
of the process ([3]) are written similarly. To perform the multi-dimensional integration we 
use the REWIAD code from the CERNLIB program library and we control the accuracy at 
the level of 2-3 percent. 



We compare t 
Collaboration 



;he results of our calculations with the experimental data from the CDF 
221 ] . The kinematic region under consideration is defined by the following 
conditions VS = 1960 GeV, 1 771,2 1 < 0.9, k 1T > 14 GeV and k 2T > 13 GeV. The results of 
our calculation are shown in the Figs. [TJI2J The description of the data can be considered 
as good at the whole kinematic region of variables M, and A(f. It is very interesting 
to compare the relative weight of different processes in connection with the a s order which 
they have. So, the contribution of RR —>■ 77 process, which is order of a 2 a 2 , is greater, by 
factor 5 — 10, than the contribution of QQ — > 77 process, which is order of a 2 . This fact 
demonstrates a dominance of the gluon induced interactions at high energies. We see that 
suppression by factor a 2 is smaller that enhancement of the gluon density over the quark 
density in a proton at the small x. The contributions of the 2 — > 3 processes (QR — ► 775 and 
QQ ~ > 119)i which are order of a 2 a s , would be unimportant for the theoretical approach 
used here. At first, the events when the diphoton and the quark (gluon) jet are produced with 
close (pseudo)rapidities are suppressed by the isolation cone condition. At second, the events 
when the diphoton and the quark (gluon) jet are produced with the large (pseudo)rapidity 
gap are included in the presented consideration via the effective vertices for Reggeized parton 
interactions in the uncollinear factorization scheme, in which we use unintegrated over the 
transverse momentum parton distribution functions. 

IV. CONCLUSIONS 

We have shown that it is possible to describe data for the prompt diphoton production 
in high-energy pp collisions at the Tevatron Collider in the leading order of the QMRK ap- 
proach. The scheme of our calculation bases on the hypothesis of quark and gluon Reggeiza- 
tion in the hard processes at high energy. Our results demonstrate that the QMRK approach 
is a powerful tool in the high-energy phenomenology. It is evidently that the QMRK ap- 



proach can be used also for description of the photon plus jet production 



20| and the photon 



plus heavy quark production 21[ at Tevatron and the LHC Collider. Our results for these 



processes will be presented in the future publications [23J 
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22]. The curve 1 is the 



FIG. 1: The diphoton px distribution from the CDF Collaboration 
contribution of the QQ — * 77 process, 2 is the contribution of the RR — > 77 process, and 3 is their 
sum. 13 




FIG. 2: The diphoton mass distribution from the CDF Collaboration 22j]. The curves are the 
same as in Fig. 1. 
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FIG. 3: The diphoton Aip distribution from the CDF Collaboration 
as in Fig. 1. 
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The curves are the same 



